Transport and Helfand moments in the Lennard-Jones fluid. I. Shear viscosity 
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We propose a new method, the Helfand-moment method, to compute the shear viscosity by equi- 
librium molecular dynamics in periodic systems. In this method, the shear viscosity is written as 
an Einstein-like relation in terms of the variance of the so-called Helfand moment. This quantity, 
is modified in order to satisfy systems with periodic boundary conditions usually considered in 
molecular dynamics. We calculate the shear viscosity in the Lennard-Jones fluid near the triple 
point thanks to this new technique. We show that the results of the Helfand-moment method are 
in excellent agreement with the results of the standard Green-Kubo method. 
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I. INTRODUCTION 

Since Maxwell's first paper [l|, 0] on the kinetic theory 
of gases, shear viscosity as well as the other transport 
properties are known to find their origin in the micro- 
scopic motion of atoms and molecules composing matter. 
However, it is only since the fifties that exact formulas are 
known to calculate the transport coefficients in terms of 
the microscopic dynamics. These so-called Green-Kubo 
formulas give each transport coefficient as the time inte- 
gral of the autocorrelation function of some specific mi- 
croscopic flux associated with the transport property of 
interest 

@, 0,1,11 • Today, the Green-Kubo technique 
allows us to calculate numerically the transport coeffi- 
cients by simulating the molecular dynamics of systems 
with a finite number of particles and periodic boundary 
conditions. 

On the other hand, Einstein classic work on Brownian 
motion showed that transport properties such as diffu- 
sion can also be understood in terms of random walks. 
It was Helfand |7j who identified in 1960 the fluctuating 
quantities which, by their random walk, are associated 
with each transport coefficients. These fluctuating quan- 
tities are the so-called Helfand moments and are the cen- 
troids of the conserved quantity which is transported. In 
principle, each transport coefficient can thus be obtained 
from the linear increase of the statistical variance of the 
corresponding Helfand moment by the so-called general- 
ized Einstein relations. Nevertheless, it is not yet known 
today how these Helfand moments should be defined in 
molecular dynamics with periodic boundary conditions, 
which limits their use in numerical simulations. 

The purpose of the present paper is to derive an ana- 
lytical expression of the Helfand moment associated with 
viscosity for molecular dynamics with periodic boundary 
conditions, and to apply the Helfand-moment method to 
the calculation of shear viscosity in the Lennard-Jones 
fluid near the triple point. 

Already in the first numerical calculation of viscosity in 
1970 Q, the algorithm of Alder et al. was based on gen- 
eralized Einstein relations derived from the Green-Kubo 



formulas. The application of the pure Green-Kubo tech- 
nique by equilibrium molecular dynamics to the Lennard- 
Jones fluid has been performed a short time afte r by 
Levesque et al and later by Schoen and Hoheisel [loj . 
However, until the middle of the eighties and the work of 
Schoen and Hoheisel [13] , as well as the one of Erpenpeck 
using the Monte-Carlo Metropolis method [ll|, nonequi- 
librium molecular dynamics was predominantly used for 
the computation of shear viscosity [l^, EH EH, [151 Il6j . 

At the end of the eighties and the beginning of the 
nineties, the generalized Einstein relations started to be 
used for calculating the transport coefficients. An alter- 
native equilibrium molecular dynamics method has been 
proposed in which the variance of the time integral of 
the microscopic flux is calculated flTl [Lsl ]. Actually it 
is the analog of the method of Alder et a/.'s Q for soft 
sphere potential systems. Recently, this technique has 
been applied by Meier et al. [l9| , and by Hess and Evans 
[2fj| . the latter having rather considered an equilibrium 
ensemble of time averages of the flux. In this context, 
two important points were discussed. The first concerns 
the so-called McQuarrie expression. In his book [2l| . 
McQuarrie presented Helfand's formula for the shear vis- 
cosity, but with a slightly different form. This difference 
implied a simplification of the expression, apparently giv- 
ing an important advantage compared to the original 
formula [22|, HU HH, The second point concerned 

the validity of the generalized Einstein relations in peri- 
odic systems [l7], EH H3, HH ■ Arguing that the periodic 
boundary conditions imply that the variance of the origi- 
nal expression of the Helfand moment is bounded in time, 
the generalized Einstein relations were considered to be 
impractical in periodic systems. In this paper, we will 
show that a generalized Einstein relation is available for 
viscosity after the addition of two terms to the original 
Helfand moment to take into account the periodicity of 
the system. The Helfand-moment method presents the 
important advantage to define shear viscosity as a non- 
negative quantity, satisfying the positivity of the entropy 
production. We have previously applied such a method 
to a system of two hard disks with periodic boundary 
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conditions [27}. In the present paper, we calculate the 
shear viscosity in a Lennard- Jones fluid near the triple 
point. We compare the results obtained by the Helfand- 
moment method with our own Green-Kubo values and 
those found in the literature. 

In addition, the Helfand-moment method plays an im- 
portant role in the escape-rate formalism. This formalism 
establishes direct relationships between the characteris- 
tic quantities of the microscopic chaos (Lyapunov expo- 
nents and fractal dimensions) and the transport coeffi- 
cients [H, [2^, [3(| UK]. A few years ago, such a relation 
has been studied for the viscosity in the two-hard-disk 
model Furthermore with the use of the Helfand 

moment, it should be possible to construct at the mi- 
croscopic level the hydrodynamic modes, which are the 
solutions of the Navier-Stokes equations. This approach 
called the hydrodynamic-mode method has been success- 
fully applied for diffusion [33|. 

The paper is organized as follows. In Section [Hi we 
outline the theoretical background of the generalized Ein- 
stein formula. Section IIIII is devoted to the presentation 
of our Helfand-moment method used for the calculation 
of the shear viscosity in this paper. In Section llVl we dis- 
cuss the so-called McQuarrie expression and the validity 
of the generalized Einstein formula for periodic systems. 
The results of the molecular dynamics simulations are 
given in Section [Vl The comparison with our Green- 
Kubo results and previous researches are done. Finally, 
conclusions are drawn in Section IVT1 



II. EINSTEIN-HELFAND FORMULA 

One century ago, Einstein theoretically established a 
relationship between the diffusion coefficient of Brownian 
motion and the random walk of the Brownian particle 
due to its collisions with the molecules of the surrounding 
fluid M : 



D = lim 

t— >oo 



[x(t) - x(0)Y 
2t 



(1) 



where x is the position of the colloidal particle and t the 
time. Thereafter, different studies led to the well-known 
Green-Kubo formula for the shear viscosity 7/ obtained 
by Green @,|3, Kubo H and Mori [|: 



lim 

jv,v-»oo k B TV 



dt^J ( ^(t)J^(0)), (2) 



where J™ is the microscopic flux associated with the 
shear viscosity r\. The viscosity coefficient is obtained 
in the thermodynamic limit where N, V — > 00 while the 
particle density n = N/V remains constant. In the fol- 
lowing, this condition is always assumed together with 
the limit AT, V — ► 00. The simplicity of Eq. |TJ) obtained 
by Einstein presents a particular interest. The extension 
of such a relation to the other transport coefficients could 



be useful. In the sixties, Helfand proposed quantities as- 
sociated with the different transport processes in order 
to establish Einstein-like relations for the transport co- 
efficients 0]. In particular, we have for shear viscosity 
that 



T) = lim — — — — 

' jv,v,t->oo 2k B TVt 



G^\t)-G^(0) 



(3) 



It can be shown [jj that the Einstein- Helfand relation is 
equivalent to the Green-Kubo formula ([2|) by defining the 
Helfand moment as: 



(4) 



if the corresponding microscopic flux is the time deriva- 
tive of the Helfand moment: 



J {ri) {t) = -G w (i). 



(5) 



We notice that the limit t — » 00 should be related to 
the thermodynamic limit N, V — ► 00. Indeed, for a fluid 
of particles confined in a finite box, the quantity (j4|) is 
bounded so that the coefficient ([3]) would vanish if the 
limit t — > 00 was taken before the thermodynamic limit 
N, V — > 00. Therefore, the number N of particles and 
the volume V should be large enough in order that the 
variance of the Helfand moment displays a linear increase 
over a sufficiently long time interval t, allowing the coef- 
ficient r\ to be well defined. The larger the system, the 
longer the time interval. It is in this sense that the limit 
N, V, t — > 00 should be considered in Eq. ([3]). 

Another remark is that, compared with the Green- 
Kubo formula ([2]), Eq. ([3j presents the advantage to 
define the shear viscosity as a positive quantity, satisfy- 
ing the conditions of non-negative entropy production. 



III. HELFAND MOMENT IN PERIODIC 
SYSTEMS 

Often, the molecular dynamics is simulated with pe- 
riodic boundary conditions. In this case, the particles 
exiting at one boundary are reinjected at the opposite 
boundary. Due to the periodicity of the system, particles 
in the simulation box may interact with image particles 
as well as the particles inside the original unit cell. As a 
consequence, the images of particle b may contribute to 
the force F ab applied by the particle b on the particle a: 



F a b — — 



du(r ab 
dr ab 



with 



Tab 



r a -r b -L0 



(a.b) 



(6) 



(7) 



where L the length of the simulation box and de- 
termines the cell translation vector [l8| . The range of the 
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interaction potential must be smaller than L/2 to guar- 
antee that the particle a interacts only with one of the im- 
ages of b in Eq. ((6| . The interacting pair is found by the 
minimum-image convention, ||r a — r b — Lfi^'^W < L/2. 
Here, we define the quantity 

U\a{t)=L0 {a - b \t) (8) 

which is the vector to be added to in order to satisfy 
the minimum-image convention. 

For a dynamics which is periodic in the box, the posi- 
tions should jump to satisfy the minimum-image conven- 
tion. As a consequence, the positions and momenta used 
to calculate the viscosity by the Green-Kubo method ac- 
tually obey the modified Newtonian equations 

S 

t = £ F <-> • <9 > 

6(#a) 

where Ari s ' is the jump of the particle a at time t s with 
HAri^l = L. We notice that the modified Newtonian 
equations ^ conserve energy, total momentum and pre- 
serve phase-space volumes (Liouville's theorem). 

Moreover, we see that the periodic boundary condi- 
tions imply that the Helfand moment of Eq. ((4]) is 
bounded and cannot be differentiated near the times t s of 
the jumps. In order to have a well-defined quantity, one 
should remove the discontinuities at the jumps, so that 
the Helfand moment can grow without bound. In order 
to do that, we add a term I(t) to the original Helfand 
moment ([JJ to get 

G<">(i)=5>«*(*)v«(t) +/(*)• (10) 

a 

According to Eq. ([5]), the time derivative of the Helfand 
moment must be the microscopic flux 

jin) {t) = £ E^l + 1 £ FmM yab , 
* — ' m 2 

defined with the position y ab of the minimum-image con- 
vention. In order to satisfy the equality |5]) in periodic 
systems, we show in Appendix [Al that the term Lit) must 
be given by 

a s 

- oE / dr F x (r ab ) L b]ay . (12) 

where both r ab and L b \ ay depend on the time r in the 
integral of the last term. We then obtain our general ex- 
pression for the Helfand moment in systems with periodic 



boundary conditions: 

N 

G^\t) = Ep~(*) »-(*) 

a=l 

a s 
1 f* 

- oE / dr F x(?ab) L b \ ay (13) 

where = p a i{t s ) is the momentum at the time of the 
jump t s and 8(t — t s ) the Heaviside step function defined 

as 

We notice that the quantity L b \ ay has discontinuous 
jumps in order to satisfy the minimum-image conven- 
tion. Let us point out that L b \ ay changes when the force 
F x (r ab ) vanishes, so that the last term varies continu- 
ously in time and does not present any jump. We notice 
thatthe last two terms of Eq. (TIB")) involves the particles 
near the boundaries of the box. The second term is due 
to the jumps of the particles to or from the neighbor- 
ing boxes. The third term concerns the pairs of particles 
interacting between neighboring cells. The Helfand mo- 
ment of Eq. (TIB"]) can be used to obtain the shear viscosity 
coefficient for systems with periodic boundaries thanks to 
the Einstein-like relation (J3J. 

IV. DISCUSSION 

Since the beginning of the nineties, some confusions 
have been propagated in the literature concerning the 
use of the mean-squared displacement equation for shear 
viscosity. First, it concerns the so-called McQuarrie ex- 
pression. On the other hand, several works have been 
done which have prematurely concluded that the mean- 
square displacement equation for shear viscosity is inap- 
plicable for systems with periodic boundary conditions. 
These confusions and criticisms are reported in partic- 
ular by Erpenbeck in Ref. [26]. Since these questions 
are central in this paper, this section is devoted to such 
problems in order to avoid any misconception. 

A. McQuarrie expression for shear viscosity 

In his well-known and remarkable book Statistical Me- 
chanics, McQuarrie [2l| reported the work achieved by 
Helfand [7|. The derivation he proposed is quite differ- 
ent but he obtained the same intermediate relation as 
Helfand, that is [35[ 

11 = Zkkvt [Xa{t) ~ ^(°)] 2 ^(*)^(°)) ■ 

(15) 
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Thereafter in his book, McQuarrie let as an exercise the 
derivation from Eq. (|15[) of the final expression which is 
printed in Ref. [21| as follows 



vttoc 2k TVt \^2^ Xa ^ Pa y^ ~ x a(°)P°v(°)Y 

^ \a— 1 

(16) 



while Helfand obtained 



?/h = hm 



1 



N,v,t^oc 2k B TVt 



N 



^ X a {t)Pay{t) ~ X a (0)p a y(0) 



0=1 



(17) 

The difference between both expressions is in the posi- 
tion of the sum over particles, and it seems that such a 
difference is due to a typing error. Nevertheless, the Mc- 
Quarrie expression (|16p at first sight presents a certain 
advantage compared to Helfand's one (|17[) . Indeed, the 
sum over the particles can come out of the average. Con- 
sequently, one would obtain a sum of averages no longer 
depending on the different particles. If Eq. (fTB"]) would 
hold, Eq. (JTHJ) could be rewritten as 



lim 



N 



N,v,t^oo 2k B TVt 



[xi(t)p ly (t) -xi(0)pi y (0)] ; 



(18) 

In other words, the McQuarrie relation seems to present 
the advantage that shear viscosity would be evaluated 
through a single-particle expression whereas Helfand ex- 
pressed the viscosity by a collective approach. 

The first time that Eq. (UHJ) has been considered was 
in the work of Chialvo and Debenedetti *22}. Without 
giving a theoretical proof of the validity of the last equa- 
tion or the equivalence with Eq. (fl7|) . they provided a 
numerical comparison between both methods and con- 
cluded that the difference between r/u and t/mq is small. 
Later, Chialvo, Cummings and Evans [H| also compared 
the McQuarrie and Helfand expressions and speculated 
an equivalencefor theoretical reasons. Thereafter, Allen, 
Brown and Masters [24J showed by comparison with their 
Green-Kubo results that the numerical calculations for 
shear viscosity obtained by Chialvo and Debenedetti [22j 
were incorrect. Moreover, Allen devoted a comment in 
Ref. (23[, and concluded that the McQuarrie expression 



is not valid and is not able to give shear viscosity 25 1. 

IS, 



This conclusion was confirmed later by Erpenbeck 
which settled the question. As aforementioned, this is 
not really surprising since Eq. (|16| seems quite clearly 
to be the result of a typing error. This discussion em- 
phasizes the fact that viscosity is a collective transport 
property, implying the intervention of all the particles. 



B. Periodic systems and Helfand- moment method 

The other point which was questioned is whether a 
mean-squared displacement equation for shear viscosity 
is useful for systems submitted to periodic boundary con- 
ditions [13, HH Ull • First, it was pointed out that the well 



known Alder et al. method initially developed for hard- 
ball systems [8j is not based on the Helfand expressions, 
but instead on the mean-square disp lacement of the time 
integral of the microscopic flux [2y] . 

The main doubt on the use of Helfand moments in 
"periodic systems comes from the fact that the original 
expression Q is bounded and would lead to a vanish- 
ing shear viscosity in the long-time limit. By this argu- 
ment, Allen concluded that the only correct way to han- 
dle G xy (t) is to write it as J Q G xy (r) dr, and express G xy 
■in pairwise, minimum-image form [l7j ]. In other words, 
the Alder et al. method would be the only valid method 
for studying viscosity, that is, through a method inter- 
mediate between the Helfand and Green-Kubo methods. 
Let us mention that this opinion was recently followed 
by Hess, Kroger and Eva ns I20I l3r| as well as by Meier, 
Laesecke and Kabelac [H, 137} having considered systems 
with soft-potential interactions. However, this does not 
preclude the possibility to modify the original expression 
(U|) of the Helfand moment in order to recover the mi- 
croscopic flux (jlip . This is precisely what we have done 
here above with our Helfand-moment method by adding 
the following two terms 

N 
a—l s 

i r* 

- / dT Fxijab) Lblay (19) 

to the original one. Albeit the first original term is 
bounded in time, the two new terms increase without 
bound in time because of the jumps and the interactions 
between the particle and the image particles (due to the 
minimum- image convention). Therefore, they can con- 
tribute to the linear growth in time of the variance of 
the Helfand moment. The Helfand-moment method we 
propose here is completely equivalent to the Green-Kubo 
formula and presents the advantage to express the trans- 
port coefficients by Einstein-like relations, directly show- 
ing their positivity. 



V. NUMERICAL RESULTS 

We carried out molecular dynamics simulations to cal- 
culate the shear viscosity by the Helfand-moment and 
the Green-Kubo methods. We use the standard 6-12 
Lennard-Jones potential 



u(r) = 4e 



(20) 



We use the reduced units defined in Table HI 

All the calculations we perform are done with the cutoff 
r c = 2.5cr. The equations of motion are integrated with 
the velocity Verlet algorithm [38| of time step At = 0.003. 
The initial positions of the atoms form a fee lattice and 
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Quantity 



temperature 
number density 
time 
distance 
shear viscosity 



Units 

' i i 



T* = 
n* = 
t* = t. 



TABLE I: Reduced units of the Lennard- Jones fluid. 



the initial velocities are given by a Maxwell-Boltzmann 
distribution. Thereafter, the system is equilibrated over 
3 x 10 5 time steps to reach thermodynamic equilibrium. 
After the equilibration stage, the production stage starts. 
At each time step, the microscopic flux (fTTj) and the 
Helfand moment (TT3"|) are calculated. Every 300 time 
units (10 5 time steps), we compute the time autocorre- 
lation function of the flux and the mean-square displace- 
ment of the Helfand moment for this piece of trajectory 
and average them with the previous results. Thanks to 
this method we can calculate with a very large statistics 
since we do not need to keep in memory the whole tra- 
jectory. Depending on the size of the systems (N = 1 08- 
1372), the number of pieces of trajectory varies between 
2000 and 6000, hence the total number of time steps is 
between 2 x 10 s and 6 x 10 8 . Statistical error is obtained 
from the mean-square deviation of the correlation func- 
tion or of the mean-square displacement on the trajectory 
pieces. 




FIG. 1: Viscosity at the phase point T* = 0.722 and n* — 
0.8442 for N = 1372. The plain line is the derivative of the 
mean-square displacement of the Helfand moment and the cir- 
cles the integral of the microscopic flux autocorrelation func- 
tion. 

We depict in Fig. [1] the time derivative of the mean- 
square displacement of the Helfand moment (fl~3|) and the 
time integral of the autocorrelation function of the micro- 
scopic flux (|lip. In Fig. [TJ the calculation is performed 



for N = 1372 atoms near the triple point at the reduced 
temperature T* = 0.722 and density n* = 0.8442. As we 
see, the two method are in perfect agreement. 
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FIG. 2: Viscosity at the phase point T* = 0.722 and n* = 
0.8442 as a function of the inverse of the number N of atoms. 
The circles are the results of the numerical simulations and 
the dashed line the linear extrapolation. 

We estimated the shear viscosity by a linear fit on the 
mean-square displacement of the Helfand moment. The 
fit is done in the region between 5 and 10 time units to 
guarantee that the linear regime is reached. We depict in 
Fig. [3 the shear viscosity versus the inverse N -1 of the 
system size. The linear extrapolation gives the following 
estimate of shear viscosity for an infinite system, 



V 



= 3.291 ±0.057 



(21) 



This result is in agreement with the previous works as 
reported in Table HT1 Indeed, the previous extrapolation 
found by Erpenbeck [U = 3.345 ± 0.068, Palmer [H 
if = 3.25 ± 0.08, and Meier et al. Ql] rj* = 3.258 ± 
0.033 are all in agreement with our result (f2"Tj) within the 
statistical error. 



VI. CONCLUSIONS 

In this paper, we propose a new method for the com- 
putation of shear viscosity by molecular dynamics. The 
Helfand- moment method is an adaptation of the Helfand 
formula ([3]) for systems with periodic boundary condi- 
tions by adding two terms (|12|) to the original expres- 
sion of the Helfand moment (j4|. The method consists 
in the calculation of the mean-square displacement of 
the Helfand moment (fTB")) . The variance of this quan- 
tity gives the shear viscosity by the generalized Einstein 
relation We have discussed its validity in the light 
of the discussions found in the literature of the begin- 
ning of the nineties. Thanks to this new method, we 
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have computed the shear viscosity in the Lennard-Jones 
fluid near the triple point. We showed that the Helfand- 
moment method gives the same results as the standard 
Green-Kubo method. Moreover, our extrapolated value 
of shear viscosity is in statistical agreement with those 
found in the literature. More than stating as an alterna- 
tive method to the standard Green-Kubo in equilibrium 
molecular dynamics, the Helfand- moment method is use- 
ful and plays a central role in the escape-rate formalism 
and the hydrodynamic-mode method. Indeed, in these 
theories, the Helfand moment allows us to put in evi- 
dence fractal structures at the micro scop ic level, which 
are related to the transport processes [28l l29l. [30L 13X1. |32| . 

We remark that the method can be similarly extended 
to the bulk viscosity £. As for the shear viscosity, two 
terms must be added to the original expression of the 
Helfand moment associated with the bulk viscosity. For- 
mally, this last coefficient is expressed as follows: 



, 4 1 

C H — v = lim 

s 3 ' N.y^oo 2k B TVt 



G (c) (t) - (G (c) (*)> 



where its Helfand moment is defined as 



(22) 



E Pax (t) X a (t) 
a 

EE^ } A ^ s) *(*-*«) 

a s 

lax 1 



(23) 



a^b ' 



and with G^(0) = 0. In the companion paper, we will 
present a similar method for the calculation of thermal 
conductivity [45| . 
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APPENDIX A: DERIVATION OF THE HELFAND 
MOMENT FOR THE SHEAR VISCOSITY IN 
PERIODIC SYSTEMS 



where we have used the modified Newton equations ([9]). 
The term implying the interparticle force F(r a fc) may be 
modified into 



^F x {r ab )y a (t) = \ ^F x {Y ab )y a + \^F x {r b a„ 



Vb- 



a^b 



a^b 



a^b 



(A2) 



Since the force F is central, wc obtain F x (r a b) 
—F x (rb a ), which implies that 



J2F x (r ab )y a (t) = J^i ;i x(r ob )(2/ -y fe ). (A3) 



a^b 



a^b 



which still differs from the corresponding term appearing 
in the microscopic flux defined with the minimum- 
image convention because y a — yb = y a b + Lb\ a y according 
to Eqs. (HI) and |JH]). Consequently, Eq. (|A1[) becomes 



dG^\t) 



dt 



a : b^a 



+ EE^w A ^ s) <^-^ 



di{t) 

dt 



(A4) 



Comparing with Eq. ([5]) , we should have 



= -EE^) A ^(^) 



~c\ E F x( V ab)Lb\ay 



(A5) 



whereupon I(t) can be expressed as 



By taking the time derivative of the Helfand moment 
(fTU)) . we have: 



dG^(t) 



dt 



E 



Pax 



EEp-w a ^ s) ^-^) 

a s 

di(t) 



+ E^^w 

a^b 



dt 



(Al) 



J (*) = -EE^ )A yi s) ^-^) 

a s 

1 f l 

- oE / dT Fxi^ab) L b \ ay ■ (A6) 
1 „ >,J„ Jo 



a,b^a ' 
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Authors 


year 


method 


' cut 




N 




An' 


Lcvcsquc et al. 1 9] 


1973 


GK (MD) 


N. C. 


0.722 


864 


4.03 


0.3 


Ashurst and Hoover [13] 


1973 


SSW 


N. C. 


0.722 


864 


3.88 


0.0026 


Al j ITT T A r\l 

Ashurst and Hoover 1401 


1975 


SSW 


N. C. 


0.722 


oo 


2.9 


0.1 


Hoover et al. [14] 


1980 


SSW 


2.5 


0.715 


864 


3.0 


0.15 




osw 




0.722 


108 


3.18 


0.1 


Levesque a ' 


1980 


GK (MD) 


N. C. 


0.728 


108 


2.97 


N. C. 










0.715 


256 


2.92 


N. C. 










0.722 


864 


3.85 


N. C. 


Pollock a) 


1980 


GK (MD) 


N. C. 


0.722 


256 


2.6 


0.1 








0.722 


500 


3.2 


0.2 


Evans [15] 


1981 


LE 


2.5 


0.722 


108-256 b) 


3.17 


0.03 


Heyes [41] 


1983 


DT 


2.5 


0.73 


500 


3.08 


0.24 


Schoen and Hoheisel [10] 1985 


GK (MD) 


2.5 


0.73 


500 


3.18 


0.15 


Erpenbeck [11] 


1988 


GK (MC) 


2.5 


0.722 


108 


2.912 


0.071 










864 


3.200 


0.160 












oo 


3.345 


0.068 


Heyes [42] 


1988 


GK (MD) 


N. C. 


0.72 


108 


3.2 


0.16 










256 


3.5 


0.18 












500 


3.4 


0.17 


Ferrario et al. [43] 


1991 


GK (MD) 


2.5 


0.725 


500 


3.02 


0.07 






0.7247 


864 


3.24 


0.10 










0.725 


864 


3.11 - 3.27 c) 


0.10 












2048 


3.24 


0.04 












4000 


3.28 


0.13 


Palmer [39] 


1994 


TCAF 


2.5 


0.722 


oo 


3.25 


0.08 


Stassen and Steele [44] 


1995 


GK (MD) 


3.4 


0.722 


256 


3.297 


N. C. 


Meier et al. [19] 


2004 


GEF 


2.5 


0.722 


108 


2.984 


0.089 








3.25 




256 


3.105 


0.093 








a n 




ou\J 


3 1 88 

O. -LOO 










5.0 




864 


3.314 


0.099 








2.5 




1372 


3.277 


0.098 








5.5 




2048 


3.224 


0.097 








5.5 




4000 


3.275 


0.098 












oo 


3.258 


0.033 


This work 


2007 


nivi 1 LVLLJ J 


9 r i 


fl 799 
u. I 


i ns 


O.UO 1 












256 


3.232 


0.060 












500 


3.270 


0.050 












864 


3.247 


0.033 












1372 


3.268 


0.055 












oo 


3.291 


0.029 



TABLE II: Results found in the literature for the shear vis- 
cosity in the Lennard- Jones fluid near the triple point. The 
reduced density equals n* = 0.8442, except for the result re- 
ported by Heyes (1988) gl] (n* = 0.848), and by Stassen 
and Steele g|] (n* = 0.8445) . Abbreviations: DT, differ- 
ence in trajectories method. GEF, generalized Einstein re- 
lation with integration of the flux. GK (MC), Green-Kubo 
results with Monte-Carlo method. GK (MD), Green-Kubo 
results with molecular dynamics. OSW, oscillatory shearing 
walls. SSW, steady shearing walls. TCAF, method based on 
the transverse-current autocorrelation functions. HM (MD) 
the present Helfand-moment method with molecular dynam- 
ics (MD). The infinite sign means that the value of the shear 
viscosity is obtained by extrapolation for iV — > oo. N. C. 
means that the value has not been communicated. 
a ' Values unpublished but communicated by Hoover et al. 

M- 

b ' No size dependence. 

c ' Values obtained for different thermostatting rates. 



